A route to robust double pulse excitability in optically injected semiconductor lasers 
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We present and analyse a three-dimensional model for a quantum dot semiconductor laser with 
optical injection. This model describes recent experimental single and double excitable intensity 
pulses, which are related to a central saddle-node homoclinic bifurcation as in the Adler equation. 
Double pulses are related to a period doubling bifurcation and occur on the same homoclinic curve as 
single pulses. The bifurcation scenario consolidating single and double excitable pulses is described 
in detail. 

PACS numbers: 42.65.Sf, 05.45.-a, 42.55.Px 



INTRODUCTION 



Semiconductor lasers with optical injection have been 
the subject of numerous studies to generate optical bista- 
bility [1], increase the modulation bandwidth of directly 
modulated lasers [2[, control the transverse mode struc- 
ture [3] or analyse the appearance of spatial solitons [J]. 
These experiments have also revealed a rich and com- 
plicated dynamical behaviour that has motivated many 

studies @,M 0, @, i, ED, EH . 

In its simplest form a semiconductor laser with optical 
injection constitutes an example of a forced non-linear 
oscillator which in addition to synchronisation can also 
display mult ist ability, exhibit instabilities, cascades of bi- 
furcations [12, [l3[ and sudden chaotic transitions [14J . It 
was also considered as a prototype for excitability since 
the phase dynamics can be described by the Adler equa- 
tion [15] which was initially introduced to model syn- 
chronisation in electrical oscillators. In this equation the 
locking/unlocking transition is associated with a saddle- 
node bifurcation on a limit cycle, and as a result two 
different types of transient can be observed depending 
on the initial condition, which in the presence of noise 
gives rise to excitable pulses. Further theoretical studies 
of semiconductor lasers with optical injection have also 
predicted the appearance of multipulse excitability asso- 
ciated with n- homoclinic bifurcations 0, [TO, [H| • 

Although excitability was experimentally observed in 
several laser systems [H, [13, IH| , it was found only re- 
cently in semiconductor lasers with optical injection [H^ . 
The experiment was carried out using a quantum dot 
laser (QDL), which displayed strongly damped relaxation 
oscillations [20[ and relatively low a-factor [2lj |. Sin- 
gle and double excitable pulses were observed for a wide 
range of parameters. 

Similar behaviour was also found theoretically in a 
four-dimensional dynamical system [19[, where the pres- 
ence of double pulses was related to a period doubling 
bifurcation 'crossing' [25[ a homoclinic curve as shown in 
Fig. [TJ This figure is a bifurcation diagram sketch which 
contains a saddle- node line SN, a period doubling curve 
PD and a homoclinic curve h. The 'homoclinic tooth', 
denned as the part of h which does not coincide with SN, 




FIG. 1: Mechanism for single and double pulse excitability. 
SN is the saddle-node line, PD is the period doubling curve 
and h is the homoclinic curve. The region LC contains stable 
limit cycle and the region LC 2 contains stable period- 2 limit 
cycle. Single and double pulses are found in the regions SP 
and DP correspondingly, the transition between these regions 
is the main point of this paper. 



plays a major role in the transition between single and 
double pulses. On the left of the tooth, the limit cycle 
LC undergoes a homoclinic bifurcation on SN, leading 
to a formation of a stable and an unstable equilibria con- 
nected by long and short heteroclines. Here, noise can 
induce the system to jump from the stable equilibrium 
to the vicinity of the unstable one, which results in a 
phase space excursion along the long heteroclinic trajec- 
tory that is shaped like a 'just broken' limit cycle. On 
the right of the tooth, the period-2 limit cycle LC 2 un- 
dergoes a homoclinic bifurcation on the same line SN, 
leading to the appearance of double excitable pulses DP 
in the presence of noise. 

This explanation differs from the one predicted earlier 
using conventional laser rate equations with optical in- 
jection [0, [l0, EH, wnere the multipulse excitability was 
predicted inside the main homoclinic tooth. The dis- 
tinguishing feature of our 'period doubling' excitability 
mechanism in comparison with the n-homoclinic scenario 
is that the regions in (Si, A) parameter space where sin- 
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gle and double pulses appear are separated by the tooth 
and are sufficiently large to make experimental observa- 
tion possible. 

Above the SN line the transition between the limit cy- 
cle LC and the period-2 limit cycle LC 2 is clearly defined 
by the period doubling bifurcation curve PD. There is, 
however, a continuous transformation between LC and 
LC 2 along the homoclinic tooth and the period doubling 
curve PD does not play a role in this transformation. 
Since a pulse resembles the shape of a 'just broken' limit 
cycle, there is no definite boundary between single and 
double pulses. 

This paper aims to describe the intricate bifurcation 
picture near the homoclinic tooth in order to understand 
the transition between single and double pulses. In the 
next section we introduce a simplified three-dimensional 
model which we use for the analysis. The third section 
describes the bifurcations of equilibria and limit cycles 
of the model as well as mentions the types of bistability, 
which were observed experimentally. Finally, section four 
focuses on the bifurcations near the homoclinic tooth, 
which are important to understand the transition from 
single to double excitable pulses. 
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FIG. 2: Qualitative sketch of (Si, A) parameter space with 
symbols indicating the stability of equilibria, where a stands 
for a stable equilibrium (attractor), r for a totally unstable 
equilibrium (repeller), and s for a saddle equilibrium. The 
shaded region is the locking region and the dark shaded region 
is the region of bistability between two locked states. 



III. EQUILIBRIA AND LIMIT CYCLES 



II. MODEL 

In order to describe the dynamical behaviour of a 
QDL with optical injection we use the following three- 
dimensional model for the dimensionless complex electric 
field E and the carrier population N (the formulation as 
well as the values of parameters are given in the Ap- 
pendix): 

E = \ ( 9 N+\E\* + 1 ~ (1 + ia)E + ^ + lAE > 

* = *{- N + P - BN N+W + l )- « 

The two main parameters are the injection strength Si 
and the frequency detuning A. We define detuning A as 
the difference between the frequencies of slave and mas- 
ter lasers. Positive detuning means that the wavelength 
of the master laser is larger than that of the free run- 
ning slave laser. Since Si and A can be easily changed 
in an experiment, these parameters constitute our bifur- 
cation parameter space. The parameters which describe 
the material properties are the differential gain factor g, 
the linewidth enhancement factor a, carrier to photon 
decay rate ratio rj, capture coefficient B, and the laser 
pump parameter is P. 

When control parameters Si and A are varied, the 
number and stability of the solutions of the dynamical 
system (pQ) generally change as well. The next section will 
discuss bifurcations of equilibria and limit cycles. The 
bifurcations presented in this paper were obtained with 
AUTO2000 [22], a software package for continuation and 
bifurcation problems in ordinary differential equations. 



The dynamical system (pQ) can have from 1 to 3 equi- 
libria depending on the region in the (Si, A) parameter 
space. Fig. [2] schematically represents equilibria bifur- 
cation curves in the (Si, A) parameter space. The thin 
curve SN is a saddle-node (also known as fold or limit 
point) bifurcation curve along which two new equilibria 
appear. The bold line HB is an (Andronov-)Hopf bifur- 
cation curve along which limit cycles are created. 

Since the system (pQ) is in the frequency frame of the 
master laser, a stable equilibrium corresponds to a locked 
operation of the laser. In the locked regime the laser 
operates at constant power and at the injected master 
light frequency. Locking can occur in regions with at least 
one stable equilibrium, i.e. in the regions ssa, sra, saa, 
and a. 

As in the experiment [19[ , two types of bistability can 
occur in the dynamical system (pQ). Bistability between 
two stable equilibria is found in the dark shaded region 
saa in Fig. El As we leave the region saa via the HB 
curve, one of the stable equilibria bifurcates into a stable 
limit cycle. Thus, the second type of bistability between a 
limit cycle and an equilibrium can be found in the region 
ssa above the black part of the HB curve. 

There are also four codimention-2 bifurcation points. 
Points C\ and C2 are cusp bifurcations, where three equi- 
libria bifurcate, and imply the presence of hysteresis (see 
section 8.2.2 in [23]). Points FH X and FH 2 are the 
codimension-2 fold-Hopf bifurcations, where the saddle- 
node and the Hopf curves are tangent (details on these 
and other bifurcations can be found in [23[). Using the 
classifications from [23[ , FH\ is a type 3 fold-Hopf bifur- 
cation and FH2 is a type 4 fold-Hopf bifurcation with 
reversed time. 



FIG. 3: The five different scenarios of disappearance of limit 
cycle LC. The torus bifurcation line TR starts at the fold- 
Hopf point FHi and connects to the period doubling curve 
PD. 

The HB curve can be divided into sub- and supercrit- 
ical parts. Above FH\ and below FH2 the Hopf bifur- 
cation curve is supercritical (shown in black in Fig. [2J 
leading to formation of a stable limit cycle. Between 
FHi and FH2 the Hopf bifurcation is subcritical (shown 
in gray in Fig. [2]) and therefore the created limit cycle 
is unstable. The stable limit cycle LC which exists for 
large values of detuning A can disappear in five different 
ways depending on the value of Si as shown in Fig. [3j 

(i) For Si above point FH\, LC disappears via a su- 
percritical Hopf bifurcation on HB. 

(ii) For Si just below point FH\, LC first loses its sta- 
bility on a torus bifurcation curve TR (which emanates 
from point FH\ ) and subsequently disappears via a sub- 
critical Hopf bifurcation on HB. 

(iii) For even smaller values of Si, LC undergoes a 
period doubling bifurcation on the PD curve before dis- 
appearing via a subcritical Hopf bifurcation on HB. 

(iv) The limit cycle LC can also disappear via a limit 
cycle fold bifurcation. 

(v) For very small values of Si the LC undergoes a cen- 
tral saddle-node homoclinic bifurcation on the SN curve, 
thus giving rise to single excitable pulses. 

Most interestingly, for the case where LC undergoes 
a period doubling bifurcation on the PD curve, the sta- 
ble newly formed period-2 limit cycle LC 2 can undergo 
a homoclinic bifurcation giving rise to double excitable 
pulses. 

Let us first notice that on the left of the tooth the orig- 
inal LC becomes homoclinic on the SN curve. However 
on the right of the tooth (where the original LC sud- 
denly creates a period-2 cycle LC 2 while passing through 
PDa curve) the LC 2 becomes homoclinic on the same 
curve SN while the original LC disappears via a sub- 
critical Hopf bifurcation on HB. On the other hand, if 
the original LC is followed in (Si, A) space very close 
along the homoclinic curve from the left side of the tooth 
to the right side, it will gradually wind up and become 



FIG. 4: Continuous transformation between the limit cycles 
LC and LC 2 along the homoclinic tooth. LC gradually winds 
up and becomes LC 2 . Limit cycles are projected onto the 
plane N = 0. 

the period-2 cycle LC 2 as shown in Fig. [4j Thus, the 
original LC and the period doubled LC 2 are in fact the 
same object. 



IV. TRANSITION FROM SINGLE TO DOUBLE 
PULSES 

This section describes the sequence of bifurcations 
which enables transition from single to double pulses. 
The main steps as we move down to the tip of the ho- 
moclinic tooth are presented in Fig. [5l The panels show 
qualitatively different 1-parameter bifurcation sketches 
versus Si, for the fixed values of A indicated in Fig. El 
On panel (a) the limit cycle is created from a homoclinic 
bifurcation at the left side of the tooth LT. Several fold 
bifurcations occur in 3b, 3a, lb, la leading to the cre- 
ation and destruction of several limit cycles. In PDa 
period doubling bifurcation occurs. The period-2 limit 
cycle disappears in 2c via a fold bifurcation with one of 
the limit cycles created in the fold bifurcation 2d. The 
second limit cycle created in the fold bifurcation 2d un- 
dergoes homoclinic bifurcation at the right side of the 
tooth RT. The period-2 limit cycle emanating from the 
period doubling bifurcation PDb disappears on the ho- 
moclinic curve h3a. In the fold point 2a two limit cycles 
are created, one of them undergoes homoclinic bifurca- 
tion in h3b while the second annihilates in the fold point 
2e with the limit cycle which appeared in the homoclinic 
bifurcation h3c. 

The overall bifurcation picture can be easily under- 
stood in terms of branches of limit cycles. Using such ter- 
minology, panel (a) has four branches. The main branch 
links HB to LT. The two period-2 branches connect 
PDa to RT and PDb to h3a. Finally, the detached 
branch binds h3b and h3c together. 

In each panel the rectangle indicates the region where 
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a reconnection will take place, leading to the next panel. 
Parts which will form a new branch in the next panel 
are shown in gray. This reconnection is associated with 
a limit cycle cusp bifurcation point which connects two 
limit cycle fold bifurcation curves. Different branches 
appear, depending on which one of the two curves par- 
ticipates in the bifurcation sketch. 

The transition from panel (a) to panel (b) is associated 
with the cusp point A in Fig. [6l and the reconnection oc- 
curs near lb. In panel (b) the main branch therefore con- 
nects HB to h3a (and stays unaffected in all subsequent 
transitions), one of the period-2 branches connects PDb 
to LT, while the two remaining branches are unchanged 
in this transition. 

Similarly, the transition from panel (b) to panel (c) 
occurs near 2c and is associated with the cusp point B in 
Fig.[6l In panel (c) one of the period-2 branches connects 
PDa to h3b and the detached branch now links RT to 
h3c. 

The transition to the next panel (d) is associated with 
the cusp point C. The parts of the branches with LT and 
h3c are exchanged so that the detached branch directly 
links LT to RT and one of the period-2 branches connects 
PDb to h3c. 




FIG. 5: Panels showing topological 1-parameter bifurcation 
diagrams of limit cycles around the tooth as A decreases. See 
Fig. [6] for the values of A corresponding to each panel as well 
as for notation. The curves represent maxima of limit cycles 
versus Si. 



The last transition to panel (e) is associated with the 
cusp point D. The link between the LT and RT becomes 
even more direct when a closed loop branch with a num- 
ber of folds separates from the detached branch in the 
course of this transition. 

The last panel (e) shows how this scenario ends. Sev- 
eral things happen here in no particular order: h3b and 
h3c join together; pairs of limit cycle folds (2b, 2a), (la, 
lc), (3a, 3c), (2d, 2e) and (4a, 4c) undergo limit cycle 
cusp or fold bifurcations and disappear; PDa and PDb 
annihilate as we move below the PD curve. Also, as 
we reach the tip of the tooth, more closed loop branches 
may detach from the branch connecting LT to RT and 
disappear. 

It is important to note that now both sides of the tooth, 
LT and i?T, are linked together and separated from any 
other structure. The sides are surrounded by infinitely 
many folds, which appear gradually on the previous pan- 
els [in. 



V. DISCUSSION AND CONCLUSION 

The transformation between single and double pulse 
excitability can be seen as the transition from the left 
side of the tooth LT to the right side of the tooth RT. 
Panel (a) in Fig. [5] shows that in order to follow from LT 
to RT at the base of the tooth, a period doubling bifurca- 
tion PDa must be encountered. As A decreases and we 
move down to the tip of the tooth (shown on subsequent 
panels), LT and RT are directly linked as they are no 
longer connected via PDa. The bifurcation analysis pre- 
sented here reveals the coexistence of single and double 
pulses on the same homoclinic bifurcation curve 'crossed' 
by a period doubling bifurcation curve. The main mech- 
anism which allows a continuous transition from period- 1 
to period-2 limit cycles (and therefore from single to dou- 
ble pulses) is a limit cycle cusp bifurcation. 
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VII. APPENDIX 

Theoretically, the dynamics of a QDL with optical in- 
jection is described by rate equations for the complex 
amplitude of electric field E, the quantum dot occupancy 
probability p, and the number of carriers per dot in the 
surrounding quantum well, N [2(J 0]. The equations 
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FIG. 6: Main bifurcation curves in (Si, A) parameter space near the homoclinic tooth. Dashed horizontal lines indicate the 
values of A for which 1-parameter bifurcation sketches were made in Fig. [5] LT and RT are the left and the right sides of the 
homoclinic tooth. h3a, h3b, h3c are the different parts of another homoclinic curve. HB is the Hopf bifurcation curve. PDa, 
PDb are the two sides of the period doubling curve PD. Number labels mark the limit cycle fold bifurcation curves and the 
letters refer to the sides of that particular fold curve. 



describing the dynamics read: 

+ iEA m +7s\/ t~~ -> 
P = -IdP + CN(1 -p)- v g a i 2 + P ~j |2 |£| 2 , (2) 
N = -~f n N+^--2CN(l-p). 

In these equations, the parameters are (typical values 
given in parentheses): 7 S (590 ns~ x ) the photon de- 
cay rate in the cavity, j n (1 ns -1 ) and 7^ (1 ns _1 ) 
the carrier decay rates in the quantum well and in the 
dot respectively, C (1.02 ps -1 ) is the capture rate from 
the quantum well into the dot, J (6.73 x 10" 10 A) the 
current per dot, a (0.6 nm 2 ) the interaction cross sec- 
tion of the carriers in the dot with the electric field, 
a (1.2) the line width enhancement factor, q the elec- 



tronic charge, the detuning A m is the angular frequency 
difference between the slave and master lasers, \E\ 2 is 
the photon density, v g (0.833 x 10 8 m/s) the group ve- 
locity, e (2 x 10 -22 m 3 ) the gain saturation coefficient, 
go (72 cm -1 ) the differential gain, uo (1.45 fs~ x ) the mas- 
ter laser frequency and S m the injected from the master 
laser energy density. 

Let us now derive a simpler model which for the given 
parameter values keeps all the essential bifurcations of 
the full QD model, in particular a closed homoclinic tooth 
with a period doubling curve 'crossing' it. In order to 
do this we first neglect the gain saturation term e\E\ 2 . 
We then note that in this regime p can be adiabatically 
eliminated as 

_ CN + v g a\E\ 2 
P ~ CN + 2v g( j\E\ 2 + ld 

The resulting reduced system is now only three- 
dimensional but, as required, preserves a homoclinic 
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tooth and a period doubling 'crossing' it: 




v g go(CN-j d ) 



2 V CN + 2\E\ 2 av a 



iA m E, 



■Id 



7s (1 + ia)E 



J 2CN(\E\ 2 cjv g + ld ) 



CN + 2\E\ 2 av a 



Id 



(3) 



Further, introducing the dimensionless variables 



F = E^2av g / ld , M = NC/-y d , h = t ls 

into Eqs. ([3j), we obtain: 
M- 1 



2 \ g M 



M 



-M + P-BM 



M 



1(1 + ia)F + y/Si + zAF, 



(4) 



where g = gov g /j s 
JC/(qj n j d ). 



In/is, B = C/7 n , P 



The Eqs. (|4|) now have unitless time ti, which is mea- 
sured in units of photon cavity lifetime 7" 1 , and the nor- 
malised electric field F, such that \F\ 2 represents the 
number of photons in the cavity in the volume 2crv g /j d . 
The two new unitless control parameters are the injec- 
tion strength Si = 2S m av g /(j d huj)^ and the detuning 
A = A m /7 S ; the latter is measured in units of inverse 
photon lifetime j s . The new constants are the differen- 
tial gain factor g = 1.01654, parameter r] = 1.695 x 10 -3 , 
parameter B = 1.02 x 10 3 , the pump parameter P = 
4.2845 x 10 3 . We now return Eqs. (|4|) to the original 
notations for the phase variables, E and TV, and obtain 
Eqs. ©. 
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